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Synopsis 

The  need  for  an  adequate  method  to  account  for  eddy  currents  in 
magnetic  materials  ha3  been  apparent  for  some  time,  especially  since  the 
development  of  high  quality  magnetic  amplifiers.  This  paper  correlates 
existing  theories  of  eddy  currents  in  saturated  iron.  In  particular  an 
approach  originally  suggested  by  A.  G.  Ganz  is  treated  in  a comprehensive 
mathematical  analysis  based  on  classical  field  theory. 

This  approach  seems  well  suited  to  account  for  the  influence  of  eddy 
currents  upon  the  transfer  characteristics  of  magnetic  amplifiers  using 
"sharp”  core  materials.  It  is  shown  that  the  method  is  useful  also  in 
more  conventional  applications,  such  as  inductive  and  conductive  heating 
of  solid  iron. 
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EDDY-CURRENT  PHENOMENA  IN  FERRO-MAGNETIC  MATERIALS 


Introduction 


Interest  in  eddy  currents  in  solid  iron  masses  has  kept  pace  with  the 
development  of  electromagnetic  devices  generally.  The  first  problems  arose  in 
the  design  of  eddy  current  brakes  far  flywheels.  Then  the  use  of  iron  wire  for 
telephone  lines  and  iron  rails  for  the  supply  of  power  to  alternating  current 
locomotives  led  to  new  problems.  It  was  recognized  that  the  saturation  of  the 
iron  was  an  important  factor,  and  many  authors  have  presented  their  theories  to 
take  saturation  into  account.  Following  the  appearance  of  Rosenberg's  work, 
others  have  made  an  academic  problem  of  substantiating  it  or  elaborating  upon  it. 
Recently  the  damping  effect  of  eddy  currents  in  the  solid  yokes  and  pole  pieces  of 
direct  current  machines  has  become  important.  Still  more  recently,  it  has  been 
recognized  that  the  saturation  effect  is  important  in  computing  the  core  losses  in 
thin  steel  or  alloy  sheets.  Representative  references  in  the  field  are  listed  in 
the  bibliography  (references  1 to  6). 


A new  method  for  computing  the  effect  of  saturation  has  been  suggested 
by  A.  G.  Ganz  (?).  The  method  seems  to  be  a significant  departure  from  other 
efforts  to  account  for  magnetic  non-linearity  in  iron.  Many  attempts  to  handle 
the  problem  mathematically  have  introduced  the  non-linearity  as  a sort  of  correction 
to  the  linear  theory,  always  keeping  in  mind  that  the  treatment  should  reduce  to 
the  linear  theory  as  a special  case.  An  almost  exact  parallel  can  be  found  in  the 
various  approaches  to  the  non-linear  problem  of  the  magnetic  amplifier.  On  the 
one  hand  there  is  the  formulation  by  fitting  some  useful  function,  such  as  the 
hyperbolic  sine  or  a finite  power  series,  to  the  magnetization  curve.  On  the  other 
hand,  a radically  different  approach  exists  in  which  all  resemblance  to  linear 
behavior  of  the  magnetic  material  is  discarded.  (These  various  analyses  of  the 
magnetic  amplifier  are  outlined  and  compared  in  reference  (8)  ). 


The  point  to  be  made  is  ttet  any  refinement  of  the  linear  theory  which 
attempts  to  account  for  non-linearity  leads  to  very  cumbersome  mathematical 
forms.  The  fresh  approach  using  what  might  be  termed  a limiting  case  for  a 
beginning,  may,  if  carefully  applied,  lead  to  a mathematical  formulation  even 
simpler  than  in  the  linear  case.  This  seems  to  be  true  of  the  method  suggested 
by  Ganz.  It  is  the  purpose  of  this  paper  to  show  a logical  transition  from  the 
linear  theory  of  eddy  currerits  in  solid  media  to  the  method  of  Ganz,  and  to  show 
the  applicability  of  the  method  in  typical  problems.  Two  other  theories,  those 
of  Rosenberg  and  of  Barth,  will  be  included  for  comparipatti 


Comparison  of  Linear  and  Limiting  Non-Linear  Theories 


The  linear  theory  will  be  reviewed  by  the  use  of  a simple  situation. 
An  infinite  half-space  with  its  surface  in  the  x,y  plane  is  excited  such  that 
the  total  magnetic  flux  carried  in  the  x direction  per  unit  of  width  in  the  y 
direction  is  Ra  35°°^ . There  is  no  y or  z component  of  flux.  The  z 

direction  extends  into  the  material  normal  to  the  x,y  plane.  The  permeability 
of  the  material  is  fx  , its  conductivity  is  cr  , ar.d  the  radian  frequency  is 
U>  « Takinglf,  the  induction,  tfO  be  the  field  valuable,  the  field  equations 
are  By  = Bz  = 0,  and  (using  MKS  units) 


= 


2>Bx 

It 
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me  solution  is 


go  £r±''l 


(2) 


where  B0,  a complex  number, 
unit  width  is 


is  the  induction  at  the  surface. 


The  total  flux  per 


(3) 


Substituting  into  (2) 


(jU  -Jjrtfr}) 


(4) 


Expressed  in  a different  form  equation  (4)  becomes 


y/iSur 

where  0mnT.  is  a real  number. 


The  behavior  of  B*  as  a function  of  time  and  of  depth  is  illustrated  in 
Figure  1.  The  upper  curve  is  the  total  flux  0(t)  while  the  lower  curves  give 
the  flux  density  as  a function  of  time  at  the  various  levels.  At  any  instant, 
the  flux  density  varies  as  a damped  3inu30id  with  depth;  thus  there  are  alternate 
positively  and  negatively  magnetized  bands  of  material.  The  existence  of  these 
alternate  bands  can  be  3een  also  by  inspecting  the  time  function  at  various  levels. 


Now  the  same  situation  is  analyzed  once  more,  but  with  the  linear 
magnetic  material  being  replaced  by  a material  having  the  non-linear  magnetic 
characteristic  shown  in  Figure  2.  This  material  is  magnetized  to  saturation  if 
the  field  intensity  is  different  from  zero,  and  it  is  possible  to  change  the  flux 
density  only  at  H - 0.  The  statement  implies  that  if  the  flux  density  is  changing, 
H must  be  zero.  However,  the  converse  statement  is  not  necessarily  true;  if  H 
is  zero,  the  flux  density  is  not  necessarily  changing.  Thus  it  is  possible  to 
have  regions  within  this  material  where  H is  zero  but  where  the  flux  density  can 
have  any  constant  value  le3s  than  or  equal  to  the  saturation  induction.  The 
particular  constant  value  would  depend  upon  the  state  in  which  the  material  was 
left  during  Some  previous  process. 


fit  <t> 


The  mechanics  of  supporting  a function  0(t) 
now  be  investigated.  First  it  will  be  assumed  that  H at  any  point  is  a periodic 
function  with  fundamental  radian  frequency  CO  and  no  average  value.  If  this 


will 


function  has  zeroes  only  at  discrete  values  of  (wi)  which  are  T radians  apart, 
the  corresponding  induction  B will  be  a square  wave  of  period  2rr/tO  and  ampli- 
tude Bs-  On  the  other  hand,  a square  wave  of  B at  any  point  can  be  supported  by 
a discontinuous  function  H ( cot),  due  to  the  possibility  (discussed  above)  that 
B may  remain  constant  when  H is  zero.  Thus  it  seem3  reasonable  to  postulate  that 
the  induction  at  any  point  in  the  material  is  either  constant  at  3ome  value  de- 
pending on  previous  treatment  or  a sauar e-wave  with  amplitude  Ba  and  period 
2ir/w. 
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Inspection  of  Figure  1 gives  seme  insight  as  to  how  the  square  waves 
giving  B as  a function  ot  time  might  combine  to  give  a total  flux  0 (t)  which  is 
sinusoidal.  In  Figure  1,  the  induction  at  any  level  s is  a sine  wave  versus 
time,  but  whose  phase  is  shifted  with  respect  to  the  level  above.  The  next 
question  is,  can  a succession  of  square  waves  which  are  phase  shifted  with  re- 
spect to  one  another  add  to  a sine  wave?  Figure  3 shows  a finite  group  of  these 
square  waves  and  their  instantaneous  sum.  (The  resemblance  to  the  calculation  of 
the  m.m.f.  of  a distributed  winding  in  a rotating  machine  is  apparent).  If 
there  are  enough  square  waves,  phase-shifted  in  a particular  manner  which  can 
be  computed,  a smooth  sine  wave  0 (t)  results. 


The  computation  of  this  phase  shift  as  a function  of  depth  is  better 
approached  by  a change  of  variable.  At  some  instant  say  ( LO  t)^  the  induction 
above  the  level  zi(Figure  3)  is  - Bs  while  below  that  level  it  is  -+*  Bs.  At  a 
later  instant  (£Ot)2  the  induction  above  the  level  Z2  has  become  -Ba  while 
below  Z2  it  remains  -f-Bs.  Thus  a surface  of  separation  between  the  two  saturated 
states  has  moved  from  Z]_  to  Z2  during  the  interval  (t2  - t]_).  The  phase  shift 
between  the  square  waves  of  induction  at  these  two  levels  will  be  oo  (t2  - t]_) 
radians.  The  distance  (z2  - that  the  surface  of  separation  has  moved,  and 
the  phase  shift  W (t£  - t]_)  radians,  are  both  related  to  the.  change  in  flux  per 
unit  width  &0  which  has  occured  during  that  time. 


The  surface  of  separation  will  stop  when  there  is  no  further  need 
to  "subtract  flux"  by  causing  this  change  to  take  place.  The  maximum  distance 
that  this  surface  travels,  beginning  at  z - 0,  will  be  the  depth  of  penetration 
S . Wien  the  surface  of  separation  has  reached  the  depth  3 « S , the  flux 
per  unit  width  will  be  - 0max.  The  depth  of  penetration  is' 


Bs 


(6) 


This  movement  from  3-0  to  S has  occurred  in  a half  cycle,  since  the 

flux  per  unit  width  nas  been  changed  from  0,^  to  -jjav  Thus,  the  phase  shift 
between  square  square  waves  of  B versus  ( w t J at  j = o an<*  3*^  is  ^ radians. 

The  d epth  of  penetration  & is  seen  to  be  a variable  depending  on 

fe,  while  in  the  linear  material  the  depth  of  penetration  is  a property  of 
material  at  a given  frequency  (equation  4).  This  behavior  is  expected, 
since  in  the  linear  material  it  is  possible  to  increase  Bx  max  if  more  flux  is 
required;  in  the  nonli.iear  material  the  induction  is  limited  to  Bs,  and  more 
flux  can  be  obtained  only  if  § increases* 


The  position  of  the  surface  of  separation  will  now  be  determined  as 
a function  of  the  flux  0.  Let  the  location  of  the  surface  be  z' . The  flux 
per  unit  width  contributed  when  the  surface  moves  from  z'  to  ( ji  4-  At,'  ), 
changing  the  state  from  *f-Bs  to  - Bs,  will  be  a ° 


-2BS  A 3' 


(7) 


! 


- u - 


Passing  to  the  limit  A(fr->0  yields  the  differential  equation  for  z' 

13'  . . .1  (s) 

The  integration  will  be  performed  beginning  at  the  time  when  (p-  + Prf\A%_ 
and  3*  ■»  o • (Calculations  begun  when  <p=  — would  proceed  in  the 

same  way  wit v a change  in  sign,  such  that  s'  is  always  a positive  in- 
creasing function.)  Thus,  j. 


’^)  =•  ZB,  )p 

*t*ax 


rW=  1 


,x 

zes  ’ 


The  next  problem  is  to  compute  the  pliase  shift  at  different  depths. 
The  square  wave  of  B versus  ( (A  t)  at  any  level  z may  be  expressed  by 


£(wi)  = Bs  □ (“t  ~ f) 


where  the  symbol  w D (**)  " represents  a periodic  function  of  x which  is  a 
square  wave  of  amplitude  1 and  period  ( 2 t/u>)  . The  angle  ^ is  the  phase 
displacement  between  the  origin  of  the  square  wave  at  depth  z and  the  origin 
of  the  square  wave  at  z = 0,  so  ^ will  be  a function  of  z,  The  origin 
of  the  square  wave  is  taken  at  the  center  of  the  positive  flat  portion. 

Equation  (10)  was  derived  assuming  that  at  time  zero,  the  flux  is 
a positive  maximum  and  is  being  decreased.  Thus,  at  time  zero  all  the  material 
was  previously  magnetized  to  -f-  Bs  induction,  and  will  be  changed  during  the 
subsequent  half  cycle  to  - Bs  induction.  At  z ■ 0 the  origin  of  the  square 
wave  would  be  located  at  — ’fl/Z.  , since  at  time  zero,  that  value  of 

( u*t  ) must  be  Ty 2,  , The  further  dependence  of  y/  upon  z can  be  de- 

duced from  Equation  (10; , since  when  the  surface  of  separation  is  located 
at  —‘J'  we  must  have  ( u>t  — y-’  ) « tr/2  at  that  level.  Thus,  setting 


i(^)  - 3- : 


ZBs 


If  now  0 = Re  is  substituted  in  (12),  the  relationship  between 

and  z for  that  case  is 

cm  (§ +y^)]  = (13) 
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Thus,  ^(0)~  -1j£  and  , signifying  that  the  square  waves  are 

phase  shifted  by  a naif  period  ever  the  depth  £ • For  , that  is, 

below  the  depth  of  penetration,  no  change  in  induction  occurs. 

The  results  of  this  computation  for  B as  a function  of  depth  and 
time  are  summarized  in  Figure  4«  The  same  values  of  and  oO  are  used 
in  Figures  1 and  4,  while  values  of  /LL  and  CT  have  been  chosen  so  that 
the  extinction  of  the  field  takes  place  in  about  the  same  depth  in  both  cases. 
The  peak  value  of  induction  is  seen  to  be  much  higher  in  Figure  1 than  in 
Figure  4,  while  the  phase  shift  is  larger  in  Figure  4*  The  different  phase 
shift  might  have  been  anticipated;  in  the  linear  case.  Figure  1,  the  contribu- 
tion to  (t)  by  the  lower  levels,  where  the  phase  is  advanced,  has  been  re- 
duced due  to  the  damping  affect.  In  the  non-linear  case  this  gradual  damping 
does  not  exist.  This  comparison  also  validates  the  idea  of  using  square  waves 
of  flux;  the  saturation  induction  is  reached  very  early  in  the  cycle  in  Figure 
1. 

An  Intermediate  Theory 


A theory  to  account  for  saturation  which  is  in  a sense  intermediate 
between  the  linear  theory  and  the  limiting  non-linear  theory  has  been  proposed 
by  Barth  (reference  3).  It  is  assumed  that  the  flux  density  varies  sinusoidally 
with  time  at  any  level,  a3  in  the  linear  theory,  but  that  the  amplitude  cf  the 
sine  wave  is  Bs  at  every  level  (except  possibly  in  that  part  of  the  material 
which  experiences  no  change  in  induction).  The  differential  equations  of  the 
field  are  used  as  in  the  linear  theory.  The  field  equations  become,  in  MKS  units. 


Curl  H “ J 


W = O 
J - 


(15) 


which  reduce  to  the  form 

d'w  & -O 
Curl  curl  H * -g 

fit 


(16) 


At  this  point  in  the  linear  theory,  the  relationship 
is  used,  to  yield  the  equations 


( d\V  B - 0 

1 — r “T1  ^ P 

divB-AB  = -pr- 

V G - 


which  can  be  reduced  to  equation  (l).  However,  in  the  present  case,  the 
transition  from  equation  (16)  to  equation  (17)  cannot  be  made.  Instead, 
it  is  found 


J 


( 


B - O 


e ■ 


cr 


’-T 

!>3 

2>t 


(18) 


- 6 - 


Before  proceeding,  it  is  necessary  to  expand  equations  (Id)  and  interpret  them 
according  to  the  assumed  field  distribution. 


Again  it  is  assumed  that  IT  has  only  and  x-coaponent  and  this  component 
varies  only  with  z._  Therefore  ^ I v-  * O is  automatically  satisfied.  It  is  now 
assumed  that  IT  and  H'are  co-linear  vectors,  although  they  bear  no  linear  re- 
lationship to  one  anot  her.  Thus  IT has  only  an  x-component.  By  symmetry,  it  is 
assumed  that  this  component  varies  only  with  z.  In  this  particular  situation, 
therefore,  div  7T®-  O also.  The  field  equations  now  reduce  to 


?1Hy_  7>8x 


(19) 


It  is  now  assumed  that  1^  and  Bx  at  any  level  z are  both  sinusoidal 
time  functions,  and  that  they  are  in  phase.  Thus 


Hx  = U ei*  e’ut 


(20) 


8*  = Bs 

where  both  and  H are  real  functions  of  z.  Substituting  into  equation  (19), 

/ . n J 


(21) 


(These  equations  neglect  the  hysteresis  effect,  while  in  Barth's  original  paper 
this  effect  is  taken  into  account  by  introducing  a constant  phase  difference 
|3  between  H and  B in  equation  (20).) 

On  performing  the  differentiation  in  equation  (21)  and  equating  real 
and  imaginary  parts,  the  following  differential  equations  are  obtained: 


H - P.fz  - 0 
2 b + H =■  n>rBs 


(22) 


where  the  dot  signifies  differentiation  with  respect  to  z.  By  inspection  of 
these  equations,  Barth  concluded  that  H must  be  parabolic  in  the  variable 
(^./q0 ) while  p must  be  logarithmic  in  the  same  variable.  The  result  ob- 
tained by  Barth  is  that 

( H=  He(i)2 


< 


f-  a Mf.) 


(23) 


3-/?  H. 


where  the  coordinate  z is  measured  from  the  level  zQ,  inside  the  material, 
positive  toward  the  interface  between  iron  and  *dr.  The  field  intensity  applied 
at  the  surface  is  H0.  The  flux  per  unit  width  turns  out  to  be 


4 (u/t  ~ H&h  XJ2) 

e J 


(24) 


from  which also. 


3. « JS 


Bs 


The  Theory  of  Rosenberg 


(25) 


Rosenberg's  treatment  (reference  2)  assumes  that  the  flux  density  is 
a sinusoidal  function,  without  phase  shift  from  layer  to  layer,  and  with  a con- 
stant amplitude  i&x  between  z 3 0 and  z = a.  The  amplitude  B,,aT  and  the 
depth  of  penetration  a are  related,  however,  the  connecting  variable  being  the 
int)6Dsity  t ho  surr^co  In  t>um j and  Bjny  are  related  through  the 

saturation  curve..  The  two  expressions  for  a corresponding  to  equations  (23)  and 
(25)  for  xs  are  (in  MKS  units; 


*-  / 


2 H* 

r B 


00 


-y 


rv\a.x 


(K  * <k*+.)c/d 


When  0 ax  is  giver,  and  it  is  required  to  find  HQ  and 

putation  is  baseaupon  the  elimination  of  a in  equations  (26): 


±, 


(26) 


the  com- 


2 H=  B 


i 


00  (T 


(27) 


The  saturation  curve  relates  Bma^  to  H0.  With  to  and  <T  fixed,  0max  be 

plotted  versus  H0.  Thus  H0  and  bmax  may  be  found,  and  either  of  equations  (26) 
determines  a. 


Rosenberg  noted  in  the  original  paper  that  the  two  most  prominent 
assumptions,  that  is,  that  the  amplitude  of  the  flux  density  is  constant  and 
that  the  phase  shift  from  layer  to  layer  is  negligible,  introduce  compensating 
errors  in  the  computation  of  eddy  current  loss,  Barth  found  that  the  inclusion 
of  phase  shift  but  the  omission  of  time  harmonies  yields  about  ths  same  results 
as  the  Rosenberg  treatment.  It  will  be  show  that  the  Ganz  treatment,  which 
takes  both  these  effects  into  account,  yields  higher  losses. 


Eddy  Currents  According  to  the  Treatment  of  Ganz 


A basic  assumption  in  G&nz's  treatment  is  that  a change  in  induction 
can  take  place  only  when  H = 0.  Therefore,  the  sum  of  all  the  field  in- 

thc  surface  of  separatc.cn  between  th^  states  —S g and  must  be 


- 8 - 


zere.  This  statement  allows  the  derivation  of  a relationship  among  the  flux,  the 
externally  applied  field,  and  the  eddy  currents. 


The  field  equations  which  relate  the  density  of  eddy  currents  to  the  rate 
of  change  of  magnetic  induction  are 


(28) 


The  second  of  equations  (28)  says  that  the  electric  field  intensity  E is  every- 
where irrotational,  because  the  induction  can  be  only  ± T B,  where  T* is  the 
unit  vector  in  the  x-direction.  The  only  possible  irrotational  vector  E under  the 
conditions  set  for  the  problem  is  a vector  varying  in  time  but  not  varying  with 
pesitien  at  a given  instant.  To  find  out  whatTf'is,  it  is  necessary  to  convert 
the  second  of  equations  (28)  to  the  integral  form,  yielding 


(j)  ~E  m A r*  =■  —A  ff 


(29) 


n j j, 


where  C boun&S  and  traverses  S in  the  usual  counter-clockwise  sense. 


The  line  integral  in  aquations  (29)  will  be  zero  if  C lies  entirely  on 
one  aide  of  the  surface  of  separation  between  -Bg  and  +Bg.  However,  if  C cuts 
this  surface,  then  the  integral  over  S7  will  have  a time  derivative  depending  upon 
the  velocity  of  the  surface  of  separation  and  the  length  of  the  segment  common 
to  bothls7 and  the  moving  surface.  Thus  a discontinuity  in  E is  expected  at  the 
moving  surface. 


The  integrals  in  equations  (29)  are  easily  evaluated  in  the  particular 
example  being  considered.  Since  B has  only  an  x-component,  the  surface  "S^is 
taken  in  the  y,z  plane.  The'  Surface  of  separation  between  states  -Bg  and  -t-Bg 
is  an  x,y  plane  moving  in  the  z-directien.  The  segment  common  to  both  these 
surfaces  will  be  Ay  meters  long.  The  electric  field  E will  have  a y-component 
only,  which  will  (at  any  instant)  take  constant  values  above  and  below  the 
surface  of  separation,  with  a discontinuity  at  the  surface.  The  situation  is 
illustrated  in  Figure  5. 


Let  the  ihduction  (in  the  x-direction)  be  -Bg  above  the  surface  and 
-+-  Bg  belew,  with  the  surface  moving  at  velocity  ( ) meters/sec.  At 

time  t,  let  IT  have  area  in  the  region  of  -Bg  Induction  and  Sj,  area  in  the 

region  of  *+•  Bg  induction.  At  time  (t  -t-  ^ t)  the  surface  vdli  have  moved  a 
distance  ( t)*(  A t)  meters,  so  that  now  there  will  bw  C 5> , -t- )(A m j} 

area  in  the  region  of  -Bg  induction  and  (S2  - area  in  the 

region  of  -pBg  induction.  The  integration  in  equations  (29)  becomes 


2BS  (d}/Jt)  (&)(&$} 

At 


(30) 
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or 


2 Bs  # 
dt 


(31) 


Equation  (31)  gives  the  discontinuity  in  the  electric  field  existing  at  the 
moving  surface  separating  -Bs  and-t-B,. 

There  can  be  no  electric  field  below  the  moving  surface,  for  if  it 
were  to  exist,  energy  would  be  required  to  sustain  the  resulting  eddy  currents, 
and  there  is  no  such  energy  available.  Therefore, 

( 6,=  *e5d|' 

( fro 

During  the  succeeding  half  cycle,  Ey  reverses  sign,  since  in  that  interval  the 
induction  is  being  changed  from  -B#  to  -l-B,  and  ( /dt  ) remains  posi- 


The  eddy  currants  flow  with  uniform  density  (at  any  instant)  in  the 
space  ( O < 3/  ) 41x1  tbere  Are  no  eddy  currents  below  the  level  z1  . The 

total  current  contained  in  the  eddy  cz  ~ ant  paths  is,  by  equations  (28)  and 
(32) 


amp/mater 


(33) 


during  the  half  cycle  in  which  the  material  is  being  changed  from  +-B,  to 
-B  . During  the  follovdng  half-cycle,  I#  changes  sign.  Since  the  field 
intensity  must  be  zero  at  the  surface  between  -Bg  and  -4-Bs,  the  applied  field 
intensity  must  be  the  opposite  of  equation  (33);  in  other  words,  the  applied 
field  must  be  - Ie  amperes  per  meter. 


An  expression  for  the  applied  field  J8  can  be  obtained  from 
equations  (33)  and  (lO)j  that  is,  for  sinusoidal  flux  variation. 


In  the  present  case,  0 - 0 y cos(<«Ot)sc  that  during  the  half  cycle  when 
0 proceeds  from  -f-  0^ax  to  “ $max 


(hT  ( \-CX\  t*)t)  Su*’  ^ 

2B5  V 


amp/meter 


(35) 
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Equation  (35  is  similar  to  Equation  (9)  of  reference  6.  The  wave  form  of  Jg 
is  repeated  during  the  following  half  cycle  but  with  reversed  sign  (it  is  con- 
venient to  establish  a new  time  origin  each  half  cycle). 

Historical  Development  of  the  Non-Linear  Treatment 


Rosenberg’s  treatment  of  eddy  current  loss  in  solid  iron  was  the  first 
attempt  of  practical  value  to  account  for  the  effect  of  saturation.  It  was 
based  on  an  hypothesis  which  fitted  quite  well  with  the  overall  test  results 
then  known.  However,  the  essential  nature  of  the  field  problem  was  not  considered, 
since  Rosenberg  neglected  both  the  phase  shift  in  the  field  from  layer  to  layer 
and  the  fact  that  the  induction  at  any  depth  must  have  prominent  time  harmonics. 
Barth  made  the  next  step  by  removing  one  of  these  assumptions;  his  work  took  the 
phase  shift  into  account  but  still  neglected  the  effect  of  time  harmonics.  As  a 
result  both  Rosenberg's  and  Barth's  treatments  yield  a sinusoidal  flux  for  a 
sinusoidal  current,  which  is  known  to  be  in  error,  especially  when  core  materials 
are  used  whose  characteristics  approach  Figure  2.  It  has  been  shown  here  that  the 
Ganz  treatment  takes  into  account  both  the  phase  shift  and  the  time  harmonics. 

The  test  results  reported  by  Hale  and  Richardson  (reference  6)  confirm  its  validity. 

Applications  of  the  Ganz  Treatment 

Applications  of  this  new  theory  have  been  restricted  so  far  to  the 
computation  of  eddy-current  loss  in  thin  sheets  (reference  6)  and  the  explanation 
of  certain  effects  in  pulse  transformer  cores  (reference  7).  To  illustrate  the 
versatility  of  the  theory,  it  will  be  used  in  two  other  problems  of  continuing 
interest.  These  are  the  computation  of  losses  in  solid  iron  conductors,  and  the 
design  of  inductive  heating  apparatus  for  use  on  a charge  of  solid  iron. 

Losses  In  Iron  Wire 


The  first  problem  is  to  deduce  the  field  configuration  within  the  wire 
under  the  assumptions  used  by  Ganz.  The  current  will  flow  in  the  axial  direction, 
and  the  wire  will  become  magnetized  in  a circumferential  sense.  Suppose  for  the 
moment  that  a steady  direct  current  flows  in  the  wire . The  current  density  will 
be  uniform,  and  application  of  Ampere's  law  shows  that  H has  a value  everywhere 
within  the  wire,  being  zero  at  the  center  only.  According  to  the  assumed  magnetic 
characteristic,  the  induction  will  be  Bs  everywhere  within  the  wire,  directed 
circumferentially.  Now  if  the  direct  current  is  reduced  to  zero  the  wire  will 
remain  in  this  magnetic  state. 

Next  suppose  that  a direct  current  of  opposite  sign  is  to  be  established. 
The  new  steady  magnetic  state  will  be  at  saturation  induction,  but  in  the  opposite 
direction.  The  change  in  induction  can  occxir  only  if  H ■ 0.  Ampere's  law  leads 
to  the  conclusion  that  there  must  be  a circular  region  within  the  wire  which  carries 
no  current,  such  that  H ■ 0 at  the  boundary  of  this  region;  the  induction  is  the  old 
value  of  (say)  - B8  in  that  region.  The  new  current  to  be  established  is  forced  to 
flow  outside  this  region.  As  the  new  current  builds  up,  the  region  of  - Bs  is  forced 
to  shrink  toward  the  center,  and  after  a time  complete  reversal  cf  the  induction 
within  the  wire  may  be  accomplished. 

Figure  6 shows  the  wire,  with  the  boundary  between  state  +Ba  and 
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- Bs  located  at  r = r' . This  moving  boundary  causes  a change  in  the  flux- 
linkage  with  that  part  of  the  wire  within  the  boundary, that  is,  r < r1, 
while  the  part  of  the  wire  in  the  region  r‘  < r < R sees  no  change  in  flux 
linkage  (except  the  change  in  the  field  outside  the  wire,  which  is  common  to 
both  regions  and  need  not  be  considered).  Since  the  axial  electric  field  is  the 
same  in  both  regions,  the  resistance  drop  in  the  region  carrying  current  can  be 
equated  to  the  drop  due  to  induced  e.a.f.  in  the  region  carrying  no  current. 


Right-handed  cylindrical  coordinates  are  used  (Figure  6).  At  time  t, 
the  surface  of  separation  is  located  at  r * r' . The  electric  field  in  volts 
per  meter,  £Z) impressed  along  the  element  whose  area  is  dS  is  to  be  computed. 
The  flux  linkage  with  that  section  at  time  t is 


= -es(r/-rs)  + 6s  (I Z-r') 


webers/aeter  (36) 


At  time  ( ■£■  «+*  A'fc  ),  the  surface  has  moved  to  (r1  -f-  A r')  and  the  flux 

linkage  has  changed  to 


(37) 


On  passing  to  the  limit,  it  is  found  that 

^=-26s  f'  (38) 

which  is  independent  of  r3.  Therefore  every  element  dS  within  the  region 
0 <|f"<  sees  the  same  electric  field  impressed,  which  is 

E3=  -28^'  (0&r  < r')  09) 

Outside  this  region  the  electric  field  Ez  is  used  up  in  ohmic  drop  only. 
(Consequently  there  is  no  reactive  voltage  due  to  partial  flux  linkages  within 
this  type  of  conductor). 


The  current  density  in  the  conducting  annulus  is 

dr' 


dt 


(r'<  rs  (?) 


The  total  current  being  carried  by  the  wire  is 

£=  tt  J(  Rz-  r'*) 


(40) 


(a) 


t 


f 

I i 

L 


r 


: t 

! ) 
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or 


Mr' 


■dt)  » -ZB.T rer(  Rz - r'2)  *£ 


*S  " w ' 0 (42) 

which  can  be  integrated  to  find  the  coordinate  r1  as  a function  of  time: 

( -2ir  8s<r  (R1-/1)  dr'  = f”  L(t)dt  (W) 

The  situation  shown  in  Figure  6 calls  for  i(t)  to  be  positive,  since  the 
induction  is  being  changed  to  +BS.  If  for  example  i(t)  is  a sinusoidal 
function,  then  equation  (43)  would  become  ( 0 $ tot  S?  v) 


§Tv83ajr/3_3Rzr/+ZR3]  = h&H-Oitoi) 


(44) 


Three  situations  can  be  recognized  in  equation  (44)*  First,  the 
maximum  current  might  not  be  sufficient  to  cause  the  surface  of  separation 
to  move  all  the  way  to  the  center.  An  equivalent  effect  is  caused  by  an 
increase  in  frequency.  Under  such  conditions,  the  depth  of  penetration  will 

be  '^R-rU 


U R 


(45) 


which  can  be  computed  from  equation  (44)  as 

jTT  Bs<r  SZ(3R-j)  = ?-I 


rw<vx 


(46) 


If  now  there  is  pronounced  skin  effect,  that  is  f?,  the  depth  of 

penetration  is 


V IT  I? 


ttR  6s  LA>(r 


(47) 


which  expressed  in  terms  of  the  maximum  apparent  surface  current  density. 
Is™*  * /Z*R)  amperes  per  meter,  is 


0 \J  oorQs 


(48) 


ootr  Bs 

The  second  situation  which  may  be  identified  in  equation  (44)  is 
that  the  current  and  frequency  might  be  such  that  the  surface  of  separation 
just  reaches  the  center  (r1  = 0)  at  the  end  of  a half  cycle.  The  current 


i 

! ■ 
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necessary  to  cause  this  condition  is  found 


from  equation  (44)* 


= 'T  ir  or 

m**  3 " 


amperes 


(49) 


The  third  situation  occurs  if  the  current  is  so  large  that  r1 
becomes  zero  before  the  end  of  the  half  cycle;  if  so,  the  wire  exhibits  a con- 
stant resistance  for  the  remainder  of  the  half  cycle,  equal  to  its  d-c  re- 
sistance. 


The  effective  resistance  may  be  computed  for  the  case  of  pronounced 
skin  effect  ( & <<  fZ.  ) . The  loss  density  as  a function  of  time  will  be 

in  general  (from  equation  41) 


-r  Z 

M . I . 


2 l2(t) 


o 

watts/meter 


O'  T Tz«r  ( R*-/1 )L 

and  the  loss  per  meter  of  length  will  be 


(50) 


pet)  = 


.1 . 

t'W 


V<r  ( R'-r'2) 


watt s/me ter 


(51) 


in  which  r*  is  a function  of  time  which  can  be  deduced  from  equation  (43). 
The  result  is 


B. 


*2(t) 


i 


^ 1 ZnRir  J iwdt ' 

If  the  current  is  a sinusoidal  function. 


•watts/meter 


(52) 


pffl s J 


Bsu> 


A iA*. 


2-rrRo"  l yr~ 


an 


watts/meter  (53) 


The  average  power  dissipated  in  heat  over  the  half  cycle  will  be 


p ^ r 

1 

4 ( 

' for 

V ZirR* 

'it 

>r 

p - p 

Bs“>  . 

4/? 

rflV  irv\^x  ^ 

2tt  Rfr 

3*rr 

2e 


dQ 


(54) 


watts/meter 


(55) 


Wien  this  power  is  expressed  in  terms  of  the  depth  of  penetration,  as 
computed  from  equation  (47),  it  becomes 


i. 


1 


3tt  ZirRio- 


watts/meter 


(56) 
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The  effective  resistance  of  the  iron  wire  is  defined  by  the  expression 

P*v  * Jrm 

- Tr  (irfe)  o:Wo9ter  (57) 

A conductor  whose  cores  section  is  not  circular,  such  as  a rail, 
can  be  treated  according  to  formula  (57)  if  the  perimeter  of  the  cross 
section  is  substituted  for  the  quantity  £■$  (^ . 

The  depth  of  penetration  (equation  48)  is  identical  with  that 
given  by  Bosenberg.  It  must  be  recognized,  however,  that  there  is  a basic 
difference  between  the  field  models  assumed  by  Rosenberg  and  by  Ganz.  This 
difference  leads  to  higher  losses  computed  according  to  the  Ganz  treatment, 
as  will  be  show  in  equation  (59) • These  similarities  and  differences  will 
now  be  examined. 


Rosenberg  assumes  that  the  current  density,  at  any  instant,  is  a 
maximum  at  the  surface  of  the  conductor  and  decreases  at  a uniform  rate  with 
depth,  becoming  zero  at  depth  S • Further,  the  current  density  at  any  level 
is  in  time  phase  with  the  current  density  everywhere  else.  Thus  the  depth 
of  penetration  and  the  depthwlse  average  current  density  are  related  according 
to 


T 


Zir  RS  J, 


rwAx 


(58) 


while  the  actual  maximum  current  density,  at  the  surface  of  the  conductor,  is 
Z Ir 


According  to  the  Ganz  treatment,  there  is  an  annular  region  at  any 
instant  in  which  the  current  density  is  uniform  vhile  the  rest  of  the  con- 
ductor carries  no  current.  The  thickness  of  this  annulus  varies  with  time 
and  the  current  density  within  the  annulus  varies  with  time  also  (equations 
43  and  40).  However,  when  the  annulus  is  thickest,  the  current  is  maximum, 
and  the  current  densitv  at  that  instant  is  the  same  as  the  average  given  by 
Rosenberg  (equation  58).  Thus  it  is  expected  that  the  Ganz  and  Rosenberg 
theories  should  give  the  same  depth  of  penetration. 

The  ratio  of  loss  per  meter  computed  by  Rosenberg  to  that  com- 
puted according  to  the  Ganz  treatment  is 


Be.n 

Z J-mAx 

/[S. 

rz  1 

J M A \ « I 

Pc  Li 

l (ZirR'He-  J' 

(Z.TT  £) 

or,  the  Ganz  treatment  yields  (4-/fr  ) = 1.273  times  as  much  loss  as 

Rosenberg’s  theory  in  the  case  of  iron  conductors  with  pronounced  skin 
effect.  The  increased  loss  is  attributed  to  the  combined  effect  of  phase 
shift  within  the  field  and  the  presence  of  time  harmonics,  both  of  which 
were  neglected  by  Rosenberg. 


/ 


Inductive  Heating  of  Iron  in  Single-Phase  Field; 


The  relationship  between  the  surface  field  intensity  and  the  total 
eddy  current  given  by  equation  (33)  for  a plane  configuration  may  be  used  for  the 
treatment  of  inductive  heating  effects  if  the  least  dimension  of  the  surface  to 
be  heated  is  much  greater  than  the  depth  of  penetration.  For  example,  if  the 
applied  field  intensity  is  a sinusoidal  function,  equation  (33)  becomes 

2BS<r  i = JsrYW  (60) 

fit 

for  each  half  cycle.  The  location  of  the  surface  of  separation  between  states 
- Bg  and  -f-B8  is 

( _ 


and  the  depth  of  penetration  is 


V __  C _ / 2 

iiM*  — 6 ~ / 5 

d \J  Bi  ^ tr 


meter 


which  is  equivalent  to  equation  (4&). 

The  current  density  in  the  conducting  layer  will  be 

T—  _ f~p  ~~T  1 dis*-  cdt  , 

J — -^/  ~ \J  tjfUJ  O'  . — ■ amp/meter* 

a \j  i~  t*jt 

and  the  ioss  per  square  meter  of  surface  area  is 

P-  1— = T^1  1 8s ^ ' f 4*2 ^ 

r 3 * ^ v t ( /rzru 

The  average  loss  is 

4ii 


Biv  ~ 


3tt 


iL  / *•>' 


watts/meter* 


(62) 


(63) 


2 

watts/meter 


(65) 


which  is  the  equivalent  of  equation  (55).  Thus  for  sinusoidal  currents,  the 
conductive  and  inductive  heating  processes  are  identical  within  the  iron,  and 
the  same  loss  formula  applies.  This  fact  has  been  noted  by  Thornton  (section 
3.1  of  reference  9). 


Practical  design  of  an  inductive  heating  installation  requires  a 
means  to  compute  the  terminal  voltage  of  the  exciting  winding.  The  real  part 
of  this  voltage  (in  phase  with  the  eliciting  current)  will  supply  the  resistive 
drop  in  the  exciting  circuit  plus  the  equivalent  re£i stance  drop  induced  by 
the  eddy  current  losses.  The  reactive  part  will  consist  of  the  reactive  drop 
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induced  by  the  eddy  currents  plus  the  reactive  drop  due  to  stray  magnetic 
.fields  (for  example,  the  field  existing  in  the  space  between  the  exciting 
winding  and  the  surface  to  be  heated).  The  real  and  reactive  parts  of  the 
voltage  induced  by  eddy  currents  will  be  computed  next.  The  other  components 
of  terminal  voltage  may  be  computed  by  well-known  methods. 


The  electric  field  induced  in  the  exciting  current  sheet  due  to 
flux  within  the  iron  is  given  by  equation  (32),  which  when  combined  with  equa- 
tion (6l)  yields 


r _ /l&JZZ'  ( (om  t>) 

V rr-  ' y 1 - <UT3  OC  t J 


(66) 


which  reduces  to 


/ 2 ftf  ( crt  ^jr) 


volts/meter 


(67) 


By  Fourier  analysis,  the  fundamental  component  is  found  to  be 

" f C?  ■ 2.  , 4 ..1 


\j 


= feSs*  i Jl  ustr  + 4 causi 


L37T 


in 


(68) 


The  phase  angle  between  the  electric  field  and  the  exciting  current  is 

_/  0.5  - 2C>.(o° 

(69) 

OnG  = 2/(5  = 0.e95 


/=  Jan'1  0*5  — 2t>.(e> 


If  there  are  N turns  in  the  exciting  coil,  each  of  length  U 
meters,  distributed  over  an  axial  length  Jt  meters,  the  formula  relating 
voltage  and  current  (exclusive  of  voltages  due  to  external  ohmic  drop  and 
stray  fields  linking  the  exciting  circuit)  is 


V=  2.  S3 


volts 


(70) 


where  V is  the  rms  of  the  fundamental  voltage,  I is  the  rms  of  the  sinusoidal 
exciting  current,  f is  the  frequency  in  cycles/sec,  and  p f/V)  is  the 
resistivity  of  the  iron  in  ohm-meters.  The  power  delivered  to  the  iron  is 


r=  2.53 


U I3/l 


H 


JejT 


watts 


(71) 


The  new  theory  provides  a means  for  computing  the  power  factor  of 
the  load  reflected  into  the  exciting  winding  by  the  eddy  currents.  Such 
a computation  was  not  possible  using  the  Bosenberg  treatment,  since  the  effect 
of  phase  shifts  within  the  field  were  not  taken  into  account . The  power 
factor  of  the  eddy-current  load  is  seen  to  be  quite  high.  The  relatively 
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low  power  factor  encountered  in  inductive  heating  installations  is  due  to  the 
reactive  voltage  of  the  unavoidably  large  stray  magnetic  fields. 

The  wave  form  of  induced  e.m.f.  as  given  by  equation  (67)  has  been 
checked  experimentally,  with  the  results  being  presented  in  Figure  7 and  Table 
1.  Sinusoidal  current  was  forced  by  connecting  an  adequate  resistance  in  series 
with  the  exciting  coil.  The  sinusoidal  current  trace  is  included  for  phase 
comparison,  and  its  calibration  varies  among  the  different  oscillograms.  The 
voltage  calibration  is  the  same  for  all. 

The  specimen  tested  was  a closed  ring  of  square  cross  section,  0.5 
inch  by  0.5  inch,  with  mean  diameter  5.375  inches,  mads  of  a particular  grade 
of  cast  iron  having  good  uniformity  and  high  carbon  content.  Accordingly, 
the  saturation  induction  is  quite  low  and  the  resistivity  is  high  being  0.68 
microhm-meter.  Magnetic  data  are  given  in  Table  2.  The  ring  was  covered 
with  two  layers  of  ordinary  black  friction  tape  and  then  closely  wound  with 
143  turns  of  AVE-10  enamel  covered  copper  wire.  A search  coil  of  12  turns 
wa3  placed  outside  the  exciting  winding.  Tests  were  made  at  60  cycles. 

The  oscillograms  show  the  characteristic  sharp  rise  in  voltage  at 
the  beginning  of  each  half  cycle,  and  gradual  decay,  as  given  by  equation  (67). 
Furthermore,  they  show  the  distinctive  behavior  expected  when  the  depth  of 
penetration  becomes  approximately  equal  to,  or  greater  than^half  the  thick- 
ness of  the  ring.  The  induced  e.m.f.  is  sharply  reduced  to  a small  value  and 
remains  there  until  the  magnetic  process  starts  again  at  the  beginning  of  the 
next  half  cycle,  the  peak  voltage  at  the  beginning  of  the  half  cycle  being 
unaffected  by  the  later  complete  magnetic  re-orientation  of  the  ring. 

The  peak  to  peak  search  coil  voltage  has  been  computed  according 
to  equation  (67),  further  modified  by  Rosenberg's  suggestion  for  determining 
Bs  as  discussed  under  "Additional  Remarks"  below.  The  depth  of  penetration 
S has  been  computed  according  to  equation  (62)  using  the  same  modification. 
These  results,  together  with  the  measured  search  coil  voltage,  are  presented 
in  Figure  8.  The  oscillograms  show  the  characteristics  of  over-excitation 
at  about  20  amperes  exciting  current,  while  the  dpeth  of  penetration  becomes 
equal  to  half  the  thickness  of  the  ring  at  24  amperes.  (It  should  be  noted 
again  here  that  all  computations  are  based  upon  the  assumptions  of  an  in- 
finite plane  configuration). 

Additional  Remarks 


A very  detailed  set  of  experiments  has  been  reported  recently  by 
Thornton  (reference  9),  in  which  the  losses  in  mild  steel  pipe  were  deter- 
mined S3  a function  of  frequency  and  of  current.  The  objective  of  the  experi- 
ments was  to  present  an  empirical  method  of  design  for  conductive  and  in- 
ductive heating  of  iron  vessels  and  pipes.  It  was  found  that  the  loss  at 
constant  frequency  varied  as  the  1.57  power  of  the  current  while  the  Ganz 
theory  predicts  an  exponent  of  1.50. 

A somewhat  empirical  way  of  determining  B , originally  proposed 
by  Rosenberg,  does  in  fact  lead  to  a loss  exponent  for  this  type  of  steel 
which  is  as  near  1.57  *3  the  accuracy  of  available  information  allows. 
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Rosenberg  proposed  that  the  value  of  Bfl  to  be  used  in  a given  problem  should 
correspond  to  the  maximum  current  carried  per  unit  periphery  of  the  conductor 
(equations  26  and  27).  If  the  magnetization  curve  for  hot-rolled  steel  sheets 
(structural  or  "mild"  steel)  as  given  in  handbooks  is  used  in  this  way  to  compute 
the  loss  exponent,  excellent  agreement  with  Thorntoris  value  is  obtained.  Reviver 
the  agreement  extends  only  over  Thornton's  range  of  test  information  which  begins 
at  about  7 oersteds,  already  well  into  the  region  where  the  induction  ceases  to 
change  rapidly  with  applied  field. 

This  empirical  method  of  assigning  Bg  may  be  justified  somewhat  by 
inspection  of  equation  (55).  l’t  is  seen  that  the  loss  varies  a3  the  square 
root  of  Ba.  The  .induction  increases  only  moderately  with  increasing  applied 
field  when  appreciable  saturation  is  present,  and  the  theory  does  not  apply 
in  the  case  of  small  saturation.  In  view  of  the  small  changes  in  Bs  to  be  ex- 
pected, and  the  further  appearance  of  Bs  under  the  radical,  the  method  for 
taking  the  variation  of  Bs  into  account  is  not  critical.  Thus  the  correlation 
between  Bs  and  field  intensity  as  obtained  from  the  saturation  curve  could  be 
expected  to  yield  reasonably  good  results.  In  substantiation,  the  mechanics 
of  the  Ganz  theory  shows  that  the  maximum  field  intensity  is  just  equal  to  the 
maximum  amperes  per  unit  of  periphery,  so  that  higher  values  of  induction  than 
that  taken  from  the  saturation  curve  would  not  be  expected. 

Conclusions 


It  has  been  shown  that  the  theory  of  alternating  fields  in  solid  iron 
proposed  by  Ganz  is  a significant  advance,  in  that  the  objections  to  Rosenberg's 
treatment  are  overcane.  Furthermore,  a logical  transistion  from  the  classical 
theory  to  the  theory  of  Ganz  has  been  demonstrated.  The  losses  computed  accord- 
ing to  the  new  theory  are  somewhat  higher  than  those  predicted  by  Rosenherg,  as 
would  be  expected  since  the  new  theory  accounts  for  both  time  harmonics  of 
current  density  and  phase  shifts  in  the  magnetic  field.  Formulae  have  been  pre- 
sented which  show  the  usefulness  of  the  new  theory  as  a design  method  for  in- 
ductive and  conductive  heating  installations. 


Table  1 


Oscillograms  of  Figure  7 


Oscillogram 

1 (upper  left) 

2 

3 

4 

5 

6 (lower  left) 

7 (upper  right) 
3 

9 

10 

11 

12  (lower  right) 


Exciting  Current,  rms 

2.5 
5 = 0 

7.5 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.0 

45.0 

50.0 


Search  coil  volts 
peak  to  peak 
0„6l 

1 AA 
x • w 

1.31 

1.58 

1.97 
2.30 
2.54 
2.76 

2.97 
3.06 
3.27 
3.30 
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Table.  2 


Magnetization  Curve  of  Cast  Iron  Ring,  by  Flux-meter 


amp/meter 

B,  weber/sq.  meter 

0 

0 

1000 

0.55 

2000 

0.76 

3000 

0.85 

4000 

0.90 

6000 

0.97 

8000 

1.01 

10000 

1.04 

12000 

1.05 

J 
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CAPTIONS  FOR  FIGURES 


Figure  1. 

Figure  2. 
Figure  3. 
Figure  4. 

Figure  5. 
Figure  6. 
Figure  7. 

Figure  8. 


Flux  per  unit  width  .and  flux  density  in  the  classical  theory  with  cot 
stant  permeability c 

Assumed  magnetic  characteristic. 

Arrangement  of  square  waves  adding  to  a smooth  periodical  function. 

Sinusoidal  flux  per  unit  width  in  saturated  solid  ion  and  square 
wave 8 of  flux  density  at  various  depths. 

Definition  of  symbols  used  in  the  plane  configuration. 

Definition  of  symbols  used  in  the  cylindrical  configuration. 

Induced  e.m.f.  of  a search  coil  wound  on  a solid  iron  ring  excited 
with  sinusoidal  current. 

Computed  and  measured  search  coil  voltage. 


5 Webers /$<?*  meter 


Figure  3 
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